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ABSTRACT 

X-ray bright quasars might be used to trace dust in the circumgalactic and intergalactic medium 
through the phenomenon of X-ray scattering, which is observed around Galactic objects whose light 
passes through a sufficient column of interstellar gas and dust. Of particular interest is the abundance 
of grey dust larger than 0.1 ^m, which is difficult to detect at other wavelengths. To calculate X-ray 
scattering from large grains, one must abandon the traditional Rayleigh-Gans approximation. The 
Mie solution for the X-ray scattering optical depth of the Universe is ~ 1%. This presents a great 
difficulty for distinguishing dust scattered photons from the point source image of Chandra , which 
is currently unsurpassed in imaging resolution. The variable nature of AGN offers a solution to this 
problem, as scattered light takes a longer path and thus experiences a time delay with respect to non- 
scattered light. If an AGN dims significantly (> 3 dex) due to a major feedback event, the Chandra 
point source image will be suppressed relative to the scattering halo, and an X-ray echo or ghost halo 
may become visible. I estimate the total number of scattering echoes visible by Chandra over the 
entire sky: N ec h ~ 10 3 (r , fb/yr _1 ), where z^fb is the characteristic frequency of feedback events capable 
of dimming an AGN quickly. 


1. INTRODUCTION 

Several lines of evidence suggest that galaxies do not 
evolve in closed boxes; a noticeable fraction of material 
can escape galaxy disks or bulges where the vast major¬ 
ity of luminous baryonic matter resides. Metals are the 
dominant tracers of star-processed material, a fraction of 
which is locked up in dust. Dust may also play an es¬ 
sential role in galaxy feedback by giving the interstellar 
medium (ISM) the opac ity necessary to launch radiation 
pressure driven winds ([Murray et al. 20051. Assessing 
the abundance of extragaiactic dust will thereby provide 
insight for the evolution of the Universe, from the pro¬ 
duction of metals by star formation, to the mechanisms 
responsible for the structure of galaxies. The purpose 
of this work is to evaluate whether the phenomenon of 
X-ray scattering can be used to trace dust in the circum¬ 
galactic (CGM) and intergalactic medium (IGM). 

Distant X-ray bright AGN proffer the opportunity to 
■robe dust in the IGM or intervening dust reservoirs 


Evans et al. 1985 Miralda-Escude 


aerels 2012p. Dust grains scatter X-ray light over ar- 


1999; Corrales & 


cminute scales, producing a diffuse halo image around 
bright point sources. The first X-ray scattering halo was 
imaged by the Einstein Observatory around GX 339- 


4 (Rolf 1983). Scattering halos are now ubiquitously 


founcfarounH objects whose light passes through a sig- 


(e.g. Smith 2008 

Valencic et al.|2009 Tiengo et al.|2010; 

Xiang et al. 201 

) and extragaiactic sources whose sight 


ISM (Vaugha n et al.||200 6). 


Historically, X-ray scattering has been treated using 
the Rayleigh-Gans approximation (RG hereafter), which 
is valid when electromagnetic waves suffer negligible 
phase shift while traversin g a d ust grain ([van de Hulst 


|1957|[Qverbeck|1965[|Hayakawa|1970||Mauche &: Goren- 


stein 1986). From the perspective of high energy light, 
the condensed matter in dust grains is roughly a ball of 
electrons, so the Drude approximation to the complex 
index of refraction can also be applied (Smith & Dwek 
1998 hereafter SD98). The RG-Drude approximation 
yield's a cross-section that is exceptionally sensitive to 
the size of the dust grains (a) and inc reasin gl y efficie nt 
for softer X -rays, ctrg oc a 4 £+ 2 (e.g. Mauche & Goren- 


stein 


1986, SD98). 


RG-Drude’s sensitivity towards grain size makes it par¬ 
ticularly tempting to use in searching for large intergalac¬ 
tic dust, for which there exists both theoretical and ob¬ 
servational motivation. Balancing the forces of radiation 
pressure, gravity, and gas drag, several works show that 
mainly dust grains >0.1 /im can reach distances a few 


rnent throuf 

h hot halo gas (Barsella et al. 1989 

Ferrara 

et al. 

*4 . 

1990 

1991 Davies et al. 

1998). Such grains are 


and altering the perceived brightness of extragaiactic ob¬ 
jects without leaving the tell-tale sign of reddening. Grey 
grains leaking into the IGM will cause systematic shifts 
for photometric measurements of Type la SNe in dark 


energy surveys, altering their figur e of merit (e.g. Virey 


& Ealetl 20071 Suzuki et al. 2012). A magnitude shift 
dm on the order of a percent is large enough to shift the 
determination of cosmological parameters by 5w ~ 5m 
and 5Qm ~ —5m 1 which is on the level of the precision 


sought after by the next generation of cosmological sur- 

veys (Zhang 2008 Menard et al.12010a 

Corasaniti 2006). 


Petrie et al. 

(12006) made the first attempt to image 


tion of the z = 4 object QSO 1508+5714. From the 
absence of a detectable halo, they concluded that the 
density of IGM dust must be < 2 x 10 -6 , relative 

to the critical density. However, as one looks towards 
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larger dust grains and softer X-ray energies where the 
RG-Drude cross-section increases dramatically, the ap¬ 
proximation begins to break down. SD98 showed that 
this occurs roughly when a Mm < E^v is no longer true. 
Cases that violate this rule-of-thumb need to use the Mie 
scat tering solution, which can be calculated nu merically 
(e.g. Bohren fc Huffman|1983 Wiscombe|1980l ). Because 
the RG-Drude approximation overestimated the scatter¬ 
ing _opacity from gr ey du st by at least a factor of 10, 


Corrales & Paerels (2012 hereafter CP12) showed that 
th e limit on the abun dance of intergalactic dust found 
by Petrie et al. ([2006) could be relaxed to < 10~ 5 . 

TKis paper follows up on the work of CP 12 by evaluat¬ 
ing scattering in the Mie regime and comparing scatter¬ 
ing halo surface brightness profiles to the Chandra point 
spread function, or PSF (Section [2]). The dust scattering 
optical depth of the universe is ~1%, which is compara¬ 
ble to the brightness of the Chandra PSF wings (outside 
of 1 — 2"). The only way to view an IGM dust scatter¬ 
ing halo is if the point source brightness is suppressed 
in some way. Since scattered light is time delayed rel¬ 
ative to the point source, AGN that suddenly decrease 
in brightness may leave behind a scattering echo. I will 
evaluate the quasar population and comment on the ob¬ 
servability of X-ray scattering echoes with respect to the 
Chandra background (Section [3]). A summary of the re¬ 
sults and review of implications, should a scattering echo 
be observed, is presented in Section [4] 


2. INTERGALACTIC SCATTERING IN THE MIE REGIME 

From Equation 19 of CP12, the optical depth of the 
Universe to X-ray scattering through a uniformly dis¬ 
tributed medium is: 

t x( z s,E) = JJ (J E ,z n d (1 + z) 2 da (1) 

where z s is the redshift of the central X-ray point source, 
E is the observed photon energy, cte j2 is the scatter¬ 
ing cross-section of a photon with energy E( 1 + z), 
H = + D m (l + z ) 3 , and rid is the comoving 

grain size distribution, presumed constant such that the 
true density at any given redshift is n d (l + z) 3 . All cal¬ 
culations in this work use standard ACDM cosmology 
with Cl m = 0.3, Da = 0.7, H 0 = 75 km/s/Mpc, and zero 
curvature. 

I derive the number density as a function of grain size 
using the constraint: 


(0.1 pm < a < 1 pm). I will investigate grain pop¬ 
ulations that are made up of purely silicate or purely 
graphite grains in order to bracket the possible solu¬ 
tions that come from a mixture of the two. In calcu¬ 
lations invoking a single grain size, the delta function 
rid oc 6(a — a g ) will be used, where a g is the grain radius 
of interest. 

Using the upper limit D)p M = 10 -5 and a point source 
at z s = 2, Figure IT] shows the Mie scattering optical 
depth of IGM dust from silicate and graphite. Di electr ic 
consta nts for both grain types were taken from Draine 
(2003 1 and used in conjunctio n with the publicly avail¬ 
able Bohren & Huffman (1983) Mie scattering code. This 
work confirms the conjecture from CP12 that the RG- 
Drude approximation is a factor of 3-10 too large when 
investigating soft X-ray scattering from grey dust, but 
is still roughly valid if the dust has a characteristic size 
around 0.1 /irn. Figure [T] also shows the Mie scatter¬ 
ing solution for a grey distribution of graphite and sili¬ 
cate dust grains, which will be used later to calculate the 
scatteringhalo intensity from intergalactic dust. 

Figure [T] demonstrates that, for the most optimistic 
conditions, only 1-5% of z = 2 quasar light will scat¬ 
ter through the diffuse intergalactic medium. This is 
roughly the same fraction of photons that will scatter off 
micro-roughness in the Chandra mirrors, which provide 
the highest imaging resolution currently available in X- 
ray instrumentation. The ability to distinguish IGM or 
CGM scattered light from the Chandra PSF and back¬ 
ground depends on a number of factors, the subject of 
the remainder of this work. To begin the investigation, 
one must calculate the intensity profile of scattered light. 

The integral for the scattering halo intensity (Ih) is 
directly proportional to the apparent flux (F a ) of the 
central point source. Following CP12, I define the nor¬ 
malized intensity d'k/dfl = Ih/F a , yielding 


d'F 

dfl 


(<*,E) 



N d (;( z ) 


(l + z) 2 

x 2 


da 


Ed 


dn 


- Ot/x 


c dz 
~H~ 


da 


( 3 ) 

where a is the observation angle with respect to the point 
source center, z s is the redshift of the source, AGis the 
grain size distribution in terms of column densiUr] with 
the dimensionless function £(z) parameterizing the dust 
spatial distribution along the line of sight, daE,z/dVl is 
the differential cross-section evaluated at scattering angle 
a/x using the energy relevant to the dust grain at z, and 
x is the dimensionless distance 


-7ra 3 p g n d da = f^ GM p c 


( 2 ) 


where p g is the grain material density and p c is the 
cosmological critical density. Typical values for p g are 
2.2 g cm -3 and 3. 8 g cm~ 3 for graphite and silicate 
grains, respectively (|Draine|2 011 ). This work also draws 
upon the canonical result of Mathis et al. (1977, here¬ 
after MRN), who found that a mix of graphite”and sil¬ 
icate grains following a power law size distribution can 
describe the extinction curves of local interstellar dust. 
Their fiducial power law, n d oc a~ 3 5 , is used through¬ 
out this work, but the end points of the distribution 
are varied to approximate an extinction curve that is 
more red (0.005 /jm < a < 0.3 pm) or more grey 


fz’ c dz/H 
fo* c dz/H 


( 4 ) 


The parameter £ = 1 in the case of uniformly distributed 
dust with constant co-moving density and £ = S(z — z g ) 
for a dust screen at redshift z g . For a full derivation of 
dust scattering in a cosmological context, see CP12. 


2.1. Scattering halo intensity from IGM dust 

The amount of dust and metals predicted from the 
star formation history of the Universe does not currently 

1 As before it is presumed that N d oc except for exami¬ 

nation of a single grain size, in which case N d oc 5(a — a g ). 
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Fig. 1.— The optical depth to X-ray scattering is plotted for a z = 2 point source and f2)p M = 10 5 . Left: The 0.5 keV Mie scattering 
solution for a population of singly sized grains is shown using graphite (dashed) and silicate ( dot-dash ) materials. Right: The scattering 
optical depth for 0.1 /im (dark curves ) and 1.0 /im (light curves) as a function of energy. The Mie solution for a grain size distribution 
with 0.1 /im < a < 1 //rn is shown for silicate (green curve) and graphite (blue curve). 


match that observed thus far. Dubbed the ‘missing met¬ 
als problem’, approximately 65% of metals in the z ~ 2 
Universe can be accounted for from observatio ns o f galax¬ 


ies an d intergalactic absorption systems (Bouche et al. 


2007). There is also a deficit between the predicted and 


observed amount of dust found in galaxy disk s a nd diffuse 


regions of galaxy halos combined (Menard & Fukugita 
2012 hereafter MF12). They predict an excess of dust 
n^ lv k 3 x 10~ 6 at z ~ 0.3, the epoch that contributes 
most to an in tergalactic scattering halo (CP12). 


Shull et al. (20141 estimate that the low redshift IGM 
contains roughly 10% of all cosmic metals, implying an 
IGM metal density ss 10~ 5 in terms of the crit¬ 

ical density. This value is in agreement with the re¬ 
sults of MF12 if the depletion factors are 30%, which 
is common. In fact, recent cosmological simulations that 
include the effect of radiation pressure driven galactic 
win ds are a ble to rep roduce observations of CGM dust 
(Menard et al. 2010b h ereafter MSFR 10) using a dust- 
to-metals ratio ot 0.29 ( Zu et al~| 2011). From here for¬ 
ward I take = 3 x 10 _e> as a fiducial mass density 

for deriving the number density of dust grains used 
in Equation [3] 

Figure [2] shows the normalized scattering halo inten¬ 
sity using z s = 2 and the Mie scattering solution for 
dctE,z/d fl in Equation w) I found very little difference 
between graphite and silicate scattering halos, so only 
graphite material is shown for scattering off uniformly 
distributed IGM dust. In the Mie regime, the singly- 
sized 1 /im dust population does not produce a signif¬ 
icantly brighter scattering halo than that caused by a 
grey dust grain size distribution. 


2.2. Scattering halo intensity from CGM dust 

CP12 evaluated various sources and reservoirs of in¬ 
tergalactic dust, including absorption line systems. This 
work considers one additional reservoir: the circumgalac- 
tic medium, which I will use to describe any dusty ma¬ 
terial residing within 300 kpc of a galaxy disk. 

There is plenty of direct observational evidence that 


dust escapes galaxy disks to pollute the CGM. Dust is 
found out to 9 kpc around starburst galax y M82, imaged 
via polarized scattering (Schmidt et al. ]T976|), polycyclic 
aromatic hydrocarbon (PAM) emissio n from a warm gas 
com ponent mixed with the superwind (Engelbracht et al. 
20061, and far-infrared emission from a cold gas compo¬ 
nent mixed in with tidally stripped material (Roussel] 


et al. 2010). Hodges-Kluck & Bregman (2014 hereafter 
1 I K I 1) also o bserve large angle U V scattering from dust 
10-20 kpc above the mid-plane of 20 early type galaxies. 

Using over 10' objects from the Sloan Digital Sky Sur¬ 
vey, MSFR10 calculate the average quasar color as a 
function of impact parameter to a foreground galaxy. 
They find a reddening trend out to distances as large as 1 
Mpc from galaxy centers, leading them to conclude that 
the mass of dust in galactic hal os is com parable to that in 
galaxy disks. Recent work by Peek et al. (20141 repeats 
the experiment with galaxy-galaxy pairs and comes to a 
similar mass estimate. 

MSFR10 calculate the average CGM extinction as a 
function of impact parameter to a nearby foreground 
galaxy. Inverting Equation 30 of their paper to solve 
for the total dust mass column yields: 


M„ « 4 X 10- 3 (jsn^) 


- 0.86 


(5) 


where r p is the impact parameter to a foreground galaxy, 
ny is the V-band extinction opacity, and h is the Hub¬ 
ble constant in units of 100 km/s/Mpc. I approximate 
the scattering from dusty CGM with an infinitely large, 
homogeneous screen at z g = 0.5 containing a dust mass 
column equivalent to an r p ~ 100 kpc sight line. Due to 
the spatial distances involved, 100 kpc is associated with 
observation angles ~ 10", which is about the visible ex¬ 
tent of intergalactic scattering halos (CP12). 

MSFR10 also found that the slope of the quasar colors 
in their study were roughly consistent with local ISM. 
However, they prefer extinction curves similar to the 
Small Magellanic Cloud (SMC) due to the fact that Mg 
II absorbers, which they believe are tracers of galaxy 
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outflow, show reddening curves that lack the 2175A fea- 
tu re often attributed to PAHs and small graphite grains 
(York et al. |2006 ). SMC extinction curves also lack the 
2175A bump. HK14 also found that the spectral energy 
distribution of UV reflection from CGM dust fits well 
with SMC-like dust in many cases. 

Following their lead, I use the simplest approxima¬ 
tion that produces an SMC-like e xtinction curve - an 
MRN distribution of silicate dust (Pei 1992) with sizes 
0.005 /im < a < 0.3 pm - to model the scattering halo 
contribution from a foreground galaxy (Figure [2]). Scat¬ 
tering off CGM dust that may be occupying the halo of 
our own Galaxy is also shown, using a screen position of 
z g = 0 and the dust mass equivalent for Ah = 10 19 cm -2 . 
Both cases are dramatically dimmer than the scatter¬ 
ing one might expect from uniformly distributed, early- 
enriched IGM, demonstrating that at least the contribu¬ 
tion from dust in the Milky Way halo should be insignif¬ 
icant. 

On the other hand, the dust distribution from a fore¬ 
ground galaxy is not an azimuthally symmetric, homoge¬ 
nous scattering screen. For quasar sight lines lying within 
1 Mpc of a foreground galaxy, MSFR.10 observe a color 
excess that spans two orders of magnitude. Since the 
scattering halo intensity is directly proportional to dust 
abundance, one might expect significant variation in the 
amplitude of the dashed red curve shown in Figure [2j 
depending on the impact parameter of the quasar sight 
line. A more detailed study of X-ray scattering from 
the CGM of a foreground galaxy requires examination 
of a non-homogeneous screen for a variety of redshifts 
and impact parameters, including the additional absorp¬ 
tion that scattered light incurs in comparison to non- 
scattered light. This is beyond the scope of this work. I 
wish only to demonstrate that X-ray scattering is more 
suited to probe a near-homogeneous, diffuse component 
of the IGM. Except for very special quasar-galaxy pairs, 
scattering from foreground galaxy dust should be less of 
a concern. 

The SMC-like dust found thus far by MSFR10 and 
HK14 produces a steep reddening signature, which is 
contrary to the hypothesis entertained by CP12 - that 
dust escaping galaxy disks needs to be grey (> 0.1 /im) 
to receive sufficient radiation pressure and survive hot 
gas. One can offer a variety of explanations for the CGM 
dust populations found thus far: (i) dust removed from 
the galaxy disk via feedback mechanisms may be subject 
to grain shattering by shocks, producing a distribution 
that is similar to that of the ISM but lacking in large 
dust grains (Jones et al. 1996); or (ii) small grains are 
trapped in galaxy halos due to an equilibrium between 
gravity and radiation pressure while grains >0.1 //in es- 


cape to greater distances (Greenberg et al. 

1987 

Ferrara 

et al. 1990 Barsella et al. 1989 Davies et al 

. 199 

5). Mote 


low-density and hot environment of the CGM can sur¬ 
vive for a Hubble time, ~ 10 10 (a/Q.l /mi) years (Draine 


& Salpeter||1979| |Draine||2011|). 
A final pomt rests in the tact 


point rests in the fact that observational bias 
of the techniques used thus far may simply prohibit our 
ability to constrain the abundance of grey CGM dust. 
Color studies relying on optical to UV wavelengths are 
relatively blind to large dust grains, and scattering cross 



Fig. 2.— (Top:) The 1 keV intensity for intergalactic scatter¬ 
ing halos around a z s = 2 source in comparison to the Chandra 
PSF template (top), and their ratios (bottom). The Mie scattering 
halo from uniformly distributed IGM dust (fi|p M = 3 X 10 -6 ) is 
shown for a grey 0.1 — 1 pm grain size distribution (thick black 
curve), which is not significantly different from the distribution 
from singly-sized 1 pm grains (thin black curve). If IGM dust 
followed a more typical SMC-like dust distribution that cuts off 
at 0.3 pm, the scattering profile brightness would be roughly an 
order of magnitude dimmer (solid red curve). A dust column con¬ 
sistent with the results of M10 ( M^ = 3 x 10 -7 g cm -2 ) around a 
z g = 0.5 foreground galaxy is simulated with a screen containing 
SMC-like dust (red dashed curve). The scattering halo associated 
with a hypothesized mass of foreground Galactic dust (Nh = 10 19 
cm -2 ), presumably from the Milky Way CGM, is also plotted us¬ 
ing an SMC-like grain size distribution (red dashed-dotted curve). 
(Bottom:) The halo-to-PSF brightness ratio for the situations de¬ 
scribed above. To obtain a robust detection of intergalactic dust, 
the PSF must be suppressed by a significant factor. 


sections are strongest for grains with a size comparabl e 
to the wavelength of interest (e.g. Andersen et al.|2013). 

Due to a lack of sufficient information, 1 will retrain 
from speculating on the mass of grey dust within a few 
100 kpc of galaxy disks. However, if the red dust found 
thus far in the CGM is indicative of that leaking into 
the IGM, Figure [2] shows that the scattering halo from 
SMC-like IGM dust is a factor of 10 dimmer than that 
produced by grey dust. 


2.3. Comparison to the Chandra PSF 

Chandra provides imaging resolution unsurpassed by 
current X-ray telescopes, confining 95% of point source 
fight to a 2 pixel (approximately 1") regionF] This res¬ 
olution is ideal for studying X-ray scattering "halos from 
IGM dust, which have half-light radii ~ 15 — 30" (CP12). 
However Figure [I] shows that we might expect only ~ 1% 
of soft X-ray light to scatter from IGM dust, making es¬ 
sential a careful understanding of the PSF wings - caused 
by scattering off micro-abrasions on the instrument mir¬ 
ror - which can easily be confused with dust scattered 

2 Chandra Proposers’ Observatory Guide (CPOG) 
http: //cxc. harvard.edu/proposer/POG/ 
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photons. To create a template for the Chandra PSF 
wings, I chose the bright Galactic X-ray binary Her X- l, 


which is used frequ ently for PSF calibration (Gaetz et al. 
2004 GaetzpOlO ). Her X-l has a very low ISM column 


(Ah ~8x lO 1 ^ cm 2 , using 
et al.|[T978j), implying T sca (l keV) 
dust. 

I fit the 0.8-1.2 keV zeroth order image from a Chan¬ 
dra HETG observation of Her X-l (Obsld 2749) with a 
model containing a beta profile for the PSF core, power 
law wings (modified with an exponential decay term at 
small radii to create a smooth core-wing transition), and 
a constant background. Figure [2] plots the background- 
free PSF model in grey. Even though the level of Galactic 
scattering expected for Her X-l is on the order of inter- 
galactic scattering for quasars, the background of Obsld 
2749 far outshines the Her X-l scattering halo and domi¬ 
nates the surface brightness profile at radii > 50". Gaetz 


Schlegel et al. 


1998 


Bohlin 


0.5% from Galactic 


(2010) use deeper grating-free observations of Her X-l to 


show that the PSF wings curve slightly at much larger 
radii, but this detail is not significant for the purposes of 
this work. Section [3] will show that the observable size of 
an IGM scattering halos is typically 10 — 20", negating 
the need to model the PSF wings any further. 

In all of the cases evaluated so far, the Chandra PSF 
outshines scattering from intergalactic dust. The bot¬ 
tom portion of Figure [2] demonstrates that the scattering 
halo surface brightness will be, at most, 20-30% of the 
PSF brightness. However, the peak in Ih/IpsF occurs 
in the area where the instrument background is more 
likely to dominate (around 100" for bright sources). The 
extremely low scattering optical depths also require un¬ 
precedented knowledge of the PSF wings, which is af¬ 
fected by a number of other systematic uncertainties, 
such as the telescope effective area and pileup. It is 
therefore highly unlikely that, in the near future, the 
Chandra PSF will be accurately calibrated to within 1% 
of the point source brightness. 

It may only be possible to view a scattering halo from 
IGM dust if the PSF is suppressed in some way. Fortu¬ 
nately, the accretion processes powering AGN are known 
to be highly variable. Since scattered light takes a longer 
path, the flux from the scattering halo is delayed by a 
characteristic timescale. Incoming scattered light is pro¬ 
portional to the average flux of the point source over a 
period of time St in the past. We might therefore see an 
X-ray scattering echo from IGM dust if an AGN has sud¬ 
denly dropped in luminosity, i.e. F ps (t) « F ps (t — St). 
An AGN that dims suddenly may also fall below an in¬ 
strumental detection threshold, leaving behind a ring of 
low surface brightness - a ghost halo - with no point 
source at the center. 

Judging from the bottom portion of Figure[2j the point 
source image might have to be reduced by a factor of 20 
to bring the scattering halo surface brightness levels on 
par with the PSF wings. To be more conservative, I 
choose a factor of 1000 decrease in AGN luminosity as a 
benchmark. In this case, scattering that would have been 
~ 1% of the PSF brightness (e.g. scattering from SMC- 
like IGM dust at 10") will instead be 10 times brighter 
than the point source image, allowing for a 3cr detection. 
This situation would allow for unambiguous confirmation 
(or ruling out) of all three of the IGM dust cases exam¬ 


ined above. The next section will evaluate the flux and 
variability requirements necessary for an AGN to pro¬ 
duce X-ray scattering echoes brighter than the typical 
Chandra background. 

3. X-RAY SCATTERING ECHOES AND “GHOST” HALOS 

Scattered light observed at angle a has an associated 
time delay, depending on the distance to the scattering 
material: 


St = 


Dr 


\l-x) 


( 6 ) 


where Dr is the radial distance to the source and x is the 
position of the scatterer al ong the sig htlin e as defined in 
Equation[4|(CP12, see al s o |Triimp er fc Schonfelde rl 1973 


MiraldaTEscude 1999). Mauche fc Gorenstein (119861 


show that the angular dependence of the KG differen- 
tial scattering cross section can be approximated with a 
Gaussian of characteristic width 


a = 10.4' 


0.1 /irn 


(AY 

V keV/ 


(7) 


Even when the a pm < -Ek e v rule-of-thumb does not 
hold, the differential Mie scattering cross section follows 
roughly the same shape of the RG approximation within 
one cr, or < 50—100" in the cases presented here (SD98). 

The majority of light scatters within 2a, so every obser¬ 
vation angle a has a scattering angle corresponding to the 
furthest distance that one can see along that sight line, 
thereby minimizing x. This condition, 0 sca = a/x < 2ct, 
yields the maximal time delay: 


sr ~ a2D R ( i a 
Otmax ~ I -L ^ 

2a 


( 8 ) 


Taking a « a for the characteristic scattering halo image 
size, the total lifetime of a 1 keV halo is on the order of 
10 4 (U/j/Gpc) years. 

3.1. Constraints on AGN flux 

Figure [2] can be used to estimate how bright a distant 
AGN must be to illuminate the intergalactic medium 
above typical background levels. For most observations, 
the Chandra ACIS background is ~ 0.1 ct s -1 per chip, 
or SHbkg ~4x 10 -7 cts s -1 arcsec -2 (ACIS instrument, 
0.5-2 keV). 1 Far enough away from the PSF core, ghost 
halos should be sufficiently bright at observation angles 
ol ss 5" with d^/dCl ~ 10 -5 arcsec -2 for IGM dust and 
~ 10 -8 arcsec -2 for CGM dust. I calculated the thresh¬ 
old flux necessary to illuminate the IGM such that the 
scattering halo in the region of 5" is brighter than the 
background, 

(j\U 

(9) 


C ps -v=r ><S£>bkg 
as 1 5" 


and capable of a 3<r detection with 50 ks of exposure time 
Tex p), 

7r(5") 2 ^ C ps T exp >9 counts (10) 

where C ps is the point source count rate in the range of 
1 keV. The resulting count rate and approximate flux 
requirements are listed in Table 1. 







































6 


TABLE 1 

AGN FLUX THRESHOLDS FOR CREATING 
OBSERVABLE SCATTERING HALOS 


Dust type 

Above SB hkg 
ct/s flux a 

3 cr detection 

ct/s flux a 

Grey IGM 
MRN CGM 

0.04 

40 

3 X 10- 13 

3 x 10“ 10 

0.2 

200 

~ io- 12 
~ 10“ 9 


a 1 keV photons were used to estimate the flux, 
which has units of erg s —1 cm -2 . 

Since the halo lifetime of 10 4 years is short compared 
to cosmological time scales, properties of the quasar pop¬ 
ulation today should describe the quasar population ca¬ 
pable of producing X-ray scattering echoes. All objects 
in the ROSAT All Sky Su rvey Bright Sources Catalog 
(RASS BSC, Voges et al. 1999), which has a limiting 
flux around 5 x 1CU 13 


erg s' 1 cm' 2 , meet the minimum 
requirement bolded in Table 1. Cross-re ferencin g RASS 
BSC with the Million Quasars Catalog (Flesch 2014), I 
found that approximately 5% of these bright quasars (248 
total) are far enough away (z > 1) to expect observable 
IGM scattering halos. I chose z > 1 because the optical 
depth of the Universe rises steeply between z = 0 and 
1, but is relatively flat for z > 1 (CP 12 ). I also exam- 
inedobjects from the rassdssagn catalog (Anderson et ah 
2007), which cross references ROSAT sources with SUSS 


quasars from DR5, providing depth rather than breadth. 
Approximately 15% of the z > 1 AGN in rassdssagn 
(156 total) are bright enough to detect IGM dust. This 
implies that, overall, there should be ~ 10 3 quasars in 
the entire sky capable of producing observable scattering 
echoes from uniformly distributed, diffuse IGM dust. 

In all the catalogs examined, no sources are sufficiently 
bright to detect CGM dust models evaluated thus far. I 
will therefore omit the CGM from further discussion in 
this work. 

Figure [3] compares the Chandra ACIS background 
to the predicted IGM scattering halo intensity using 
the minimum threshold from Table 1 and a value ten 
times brighter than the minimum RASS BSC flux: 5 x 
10“ 12 erg s _1 cm -2 . There are seven z > 1 RASS BSC 
quasars that meet this higher threshold. The bottom 
portion of Figure [3] shows the maximum time delay as a 
function of observation angle and grain size. The observ¬ 
able scattering echo, which is confined to a < 5 — 20", 
should last 1 — 100 years. 

3.2. Constraints on AGN variability 

Hour to day-long activity, such as an X-ray afterglow 
to a gamma ray burst, produces an echo that probes 
only a tiny fraction of the total dust column. Therefore, 
burst activity will not make a particularly good probj^of 
the d iffuse IGM or intervening galaxies (Miralda-Escude 
1999). To produce observable scattering echoes, the cen¬ 
tral X-ray point source needs to undergo a period of sus¬ 
tained brightness before dimming by a factor > 1000. 

Simulations of AGN accretion and feedback show that 
the supermassive black holes (SMBHs) powering quasars 
rarely accumulate mass at a steady rate. They may 
spend the majority of their time in a sub-Eddington 
phase, with near-Eddington a ccreti on events t riggered by 
the sudden infall of new gas (Gabor & Bournaud 2013). 
In simulations, periods of sustained activity are punctu¬ 



Fig. 3.— Top: Point source subtracted scattering halos from 
a uniform IGM, using the same cosmological parameters as Fig¬ 
ure]^ The halo brightness is directly proportional to the average 
AGN flux (labeled F src in the plot). Scattering halos from two flux 
thresholds (black and grey curves ) are plotted in comparison to the 
typical Chandra detector background (dotted line). Bottom : The 
maximal time delay as a function of observation angle, using two 
limiting grain sizes. Scattering echoes visible by Chandra will last 
~ 1 — 100 years. 


ated by major feedback events that cau se the A GN to 


dip in luminosity by factors of 10 2 — 10 6 (Hopkins et al. 


2005 Novak et al. 2011). 


Motivated by the simulations, rapid dimming after a 
prolonged period of AGN activity is most likely to occur 
via some feedback mechanism, though technically any 
burst of activity longer than 10 3 years will do. I will call 
z/fb the characteristic frequency for such a cataclysmic 
event to occur. A rough estimate for the number of scat¬ 
tering echoes in the sky is therefore 


N e ch ~ St max I'fb N q (F > F t h,z > 1) (11) 

where N q is the total number of high- 2 : quasars with 
an apparent flux brighter than the threshold flux Eth- 
The above equation requires that the characteristic peri¬ 
odicity of the AGN growth and feedback process does 
not occur on a shorter timescale than the echo, i.e. 
itmax < zqjj 1 • Otherwise, the point source image will not 
remain suppressed relative to the scattering image. This 
work also contains the assumption that AGN brighten 
slowly so that they remain at a low luminosity for a time 
long er th an 5t max , while feedback events are rapid (Sec¬ 
tion [O]). 

The number of quasars able to illuminate the IGM, 
N q (F > F t h) is strongly dependent on a number of sit¬ 
uational factors, including the instrument background 
and image exposure time. In addition, St max and N q are 
linked because the brighter quasars will produce the most 
visibly extended halos. For the sake of theory, I examine 
the abundance of scattering echoes one might expect for 
a perfect instrument with zero background. Since the 
inner part of the halo will vanish first, a scattering echo 
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Fig. 4.— The fraction of AGN (N ec h/N q ) with X-ray scattering 
echoes of a given size, for a particular v ft,. The plotted contours, 
from black to lightest grey, are 100%, 50%, 10%, and 1%. Beyond 
the 100% contour, X-ray scattering echoes will appear as nested 
rings. 


will qualitatively appear as a ring of surface brightness 
with an angular radius corresponding to the time delay 
dim ax after the feedback event. Figure [4] plots contours 
of iV ec h/-/V q using z/fb and the halo radius calculated from 
dimax as seen in the 0.1 (im curve in the lower portion 
of Figure [3] Given a characteristic feedback timescale, 
the contours show the fraction of AGN that will have 
ring-like echoes of a particular radius. For example, if 
^fb = 10 -4 yrs -1 , then 10% of AGN will have echoes 50" 
in radius, and 50% of them will have larger but dimmer 
echoes about 100" in radius. The 100% curve in Figure [4] 
marks the point at which dt max = t/jj 1 . Larger values 
will create a series of nested X-ray echo rings, analogous 


to th ose seen around variable objects (e.g. Tiengo et al. 
2010 ). 


"From the catalog searches described in Section [3T 
N q (F > 3 x 10 -13 ,z > 1) ~ 10 3 

and 

N q (F > 5 x 10 -12 , z > 1) ~ 7. 

With df max values of 1 year and 100 years for the two flux 
thresholds, respectively, I take N q 5t max to be roughly 
constant at 10 3 years. To achieve N ec h > 1, the char¬ 
acteristic feedback frequency must be at least one major 
event per 1000 years, or vr, > 3 x 10 -11 Hz. 

The characteristic frequency of interest to this work 
skirts the limit of human timescales. Long-term X-ray 
observational campaigns covering periods up to several 
years (v > 10~ 9 Hz) revea l variations on the order of 


10 - 500% ( 

McHardy 

2013 

Soldi et al. 

2014 

Lanzuisi 

et al. 2014 

)• 

in these studies, AGN flicker, but they rarely 


burn out. Except for short-lived flares, 3 dex changes in 
luminosity are not typically observed. 

There are, however, theoretical explanations to ob¬ 
served AGN phe nomena that give hints to larger am¬ 


plitude variation. Hickox et al. (2014) produce a proba¬ 
bility distribution over time for Lagn between 10 -6 and 
one L^dd that can explain the relatively weak correlation 
between AGN luminosity and the star formation rate of 
the host galaxy as a function of redshift. Intermittent 
AGN jet activity, occurring every 10 3 — 10 5 years, may 
also explain extended radio galaxy morphologies and the 


appearanc e of young radio sources embedded in older ra¬ 
dio r elics (Reynolds & Begelman 1997 Owsianik et al. 
19981. Similar outbursts can be triggered by ionization 


instabilities in the accretion disk, which prod uce 3-4 dex 
lumin osity swings on 10 6 year timescales (Siemiginowska 
et al. 1996). Radiation pressure instabilities, which can 
describe the activity observed in Galactic microquasars, 
would lead to AGN outbursts at 10 2 — 10 5 year frequen- 
cies, depen ding o n the black hole mass and accretion rate 
(Siemiginowska et al. 2010). The above results imply a 
range of values that yield A^ ech <10. 

3.3. Constraints from accretion disk physics 

The discovery of an X-ray scattering echo is con¬ 
strained by yet another piece of physics: how quickly 
can an AGN turn off? For an X-ray scattering echo to 
be visible, the source of X-rays must dim on a timescale 
that is shorter than <5f max (< 100 years). 


Linder the clas sic thin disk accretion model (Shakura 
& Sunyaev 1973), the viscous timescale f v isc i 

describes how"quickly material at radius R will fall onto 
the central object, i.e. how long it will take the accretion 
disk to clear out. The sound speed c s and scale height of 
the disk h provide upper limits to the viscosity, v v i sc ~ 
ac s h , with a < 1 being a catch-all parameter describing 
the contributions of magnetic fields and turbulence to 
the transport of angular momentum in the disk. Using 
some values typical of supermassive black holes, 


tv 


RV 2 Me T r' 


( 12 ) 


where Mg = M/10 6 M 0 , Tg = T/10 6 K, and R s is the 
radius in units of the Schwarszchild radius. For reference, 
the radius associated with the dynamical timescale on the 
order of 5t max is 10 5 Schwarzschild radii for a 10 6 M 0 
black hole, yielding t v i sc ~ 800 years. 

To satisfy the relation f v ; sc < <5f max , a significant por¬ 
tion of the disk must be removed, leaving only material 
within R s < 40 M 6 " 1 (<5t max /100 yrs). This inequality 
also implies that black holes larger than 10' M 0 will 
not fade quickly enough to produce visible X-ray scat¬ 
tering echoes, at least under the alpha prescription for 
thin accretion disks. It is more likely that disk instabil¬ 
ities, which propagate on timescales much smaller than 
the steady-state viscous timescale, will trigger a change 
fast enough to leave behind a visible X-ray scattering 
echo. The discovery of a ghost halo from the IGM would 
thereby place interesting limits on accretion disk physics 
and the nature of feedback events responsible for turning 
an AGN “off.” 

There is additional observational evidence that AGN 
can switch from a high- to low-luminosity state in a short 
period of time. A UV light echo in the vicinity of galaxy 
IC 2497 (“Hanny’s Voorwerp”) contains fossil evidence 
of a nearby quasar that has dimmed by o ver 2-4 orders 
of magnitude within the last 70,000 years (Lintott et al. 


2009 


Schawinski et al. 


2010). Observations of Fe K fiu- 


orescence near the Galactic Center, likely arising from 
a flare up of Sgr A*, indicate that our galaxy’s central 


100 years (Ponti et al. 

2010 

Capelli et al. 

2012 

Gando 

Ryu et al.|2U12). Galactic microquasars, powerec 

by bi- 

nary accretion, can change between soft and hard states 

























































within a matter o f hours (e.g. Greiner et al.|[T996 Grin- 
berg et al.||2013). While these changes do not involve 
tactors of 1000 m luminosity, they are indicative of the 
variability achievable by accreting black holes. If these 
timescales are roughly proportional to black hole mass, a 
10 6 —10 7 M 0 AGN might change states in 10 - 100 years. 

4. CONCLUSION 

Implementing Mie scattering for a Universe suffused 
with « 3 x 10 -6 in grey dust, I have shown that 

only ~ 1% of X-ray light from a high-z source will scat¬ 
ter, not 10 — 30% as the RG-Drude approximation sug¬ 
gests. I also found that X-ray scattering off of SMC- 
like dust in the CGM of z ~ 0.5 foreground galaxies 
is much less likely to be observed in comparison to a 
dusty IGM. That’s because uniformly distributed dust 
produces more centrally concentrated scattering halos. 
However, the surface brightness from IGM dust scatter¬ 
ing would still be dimmer than the Chandra PSF wings, 
making it very difficult to disentangle an intergalactic 
scattering halo from light scattered off micro-roughness 
in the Chandra mirrors. 

A natural way to get around this problem is to search 
for a scattering echo from an AGN that has dimmed sig¬ 
nificantly, by a factor of 1000 or more. In this case, 
the point source image will be suppressed relative to the 
brightness of the scattering echo, which is proportional 
to the time-averaged flux of the AGN over <5f max (Equa¬ 
tion [sj) . An AGN that has effectively turned “off” will 
not appear as a point source, and a bare scattering echo, 
or ghost halo, will be left behind. Ghost halos will have 
observable lifetimes ~ 1 — 100 years, relative to the typ¬ 
ical soft X-ray background on Chandra. 

A catalog search of X-ray bright AGN found hundreds 
of sources sufficiently luminous to produce 3 — 10" ghost 
halos that are brighter than the Chandra ACIS back¬ 
ground. These results imply a total number of scattering 
echoes iV ec h ~ 10 3 (r , tb/yr _1 ) over the entire sky, where 
the frequency of major feedback events z'fb is by far the 
greatest source of uncertainty. 

To leave behind a scattering echo observable by Chan¬ 
dra , the contributing feedback event must meet sev¬ 
eral criteria: (i) prolonged accretion for longer than 
ifmax that produces an average apparent flux > 3 x 
10” 13 erg s -1 cm -2 , (ii) a sudden drop in luminosity 
by a factor > 1000, and (iii) the time it takes the AGN 
to dim must be shorter than St m ax . The second criterion, 
dimming by a factor of 1000, is also necessary for an un¬ 
ambiguous confirmation or ruling out of the hypothesized 
dust density nf M = 3 x 10 6 . 

Simulations of AGN accretion and feedback indicate 
that conditions (i) and (ii) can be satisfied easily, but 
a rough c ount of eligible events yields t/fh ~ 10~ 9 — 


10 7 yr 1 ( Hopk ins et al.|2005 


g 


Novak et al. 2011 Gabor 


& Bournaud|2UI3). This would imply A ec h << I. Unthe 
other hand, Intermittent jet activity invoked to explain 
some observational properties of AGN host galaxies and 
radio lobes could satisfy conditions (ii) and (iii). Sev¬ 
eral disk instability mechani sms pro posed t hus far imply 
vn, ~ 10~ 6 — 10~ 2 yr -1 (Reynolds fe Begelman 1997 


Siemiginowska et al. 1996 2010|). This produces some 
hope tor A ec h ~ 1. 

The events forming a ghost scattering halo may be in¬ 
credibly rare, but the discovery of such an object would 


have interesting implications for the fields of AGN vari¬ 
ability, feedback, and galaxy evolution. One of two 
conclusions may be reached if 7V ec h > 1: first, that 
^ GM is larger than previous expected, thereby increas¬ 
ing N q (F > F th ); or second, that the frequency of AGN 
flickering between high- and low-luminosity accretion 
states is much larger than expected (z/fb > 10 -2 yr -1 ). 
The latter conclusion would also support theories of ac¬ 
cretion beyond the thin alpha-disk prescription and help 
constrain the type of disk instabilities responsible for jets 
and other flaring activity. 

The smoking gun for a ghost halo - a diffuse, hollowed- 
out ring of soft X-ray light about 10" across - could be 
imaged with Chandra, which has sub-arcsecond imaging 
resolution. The XMM EPIC detectors have a slig htly 


lower back ground level than Chandra ACIS (Katayama 
et al.|2004 ) but a 6" FWHM, which might only barely re¬ 


solve a ring structure. All other observatories with com¬ 
parable sensitivity to soft X-rays, e.g. Suzaku or ASCA, 
lack the spatial resolution necessary to identify a ring. A 
ghost halo imaged with these instruments might resem¬ 
ble a very soft point source, with a spectrum that falls 
swiftly with increasing energy, i.e. a photon index T > 2. 

Future generations of X-ray observatories are likely 
to achieve background r ates ab out 5-10 times lower in 
comparison to Chandra (Lotti et al. 2014), decreasing 
the flux threshold and thereby improving the likelihood 
of finding an X-ray scattering echo. One of the primary 
science goals of the Athena-h mission is to detect 
AGN in the z > 6 Universe, in order to understand 
format ion mechanisms of th e earliest supermassive black 
holes (Nandra et al. 2013). Massive 10 9 M 0 SMBHs 
discovered thus far at z > 6 indicate that the earliest 
black hole seeds must have either grown at an extremely 
rapid pace or collapsed spontaneously. Either way, the 
mechanisms responsible for forming the most massive 
high redshift SMBHs are not similar to the accretion we 


see from nearb y AGN of moderate masss (e.g. Tanaka & 


Haiman 2009). Along with the instrumental increase m 
hux sensitivity and 5" imaging resolution, the increased 
IGM column of the high-z Universe might allow X-ray 
scattering echoes to contribute to the Athena+ science 
goals by probing the variable nature of young SMBHs. 
Finally, the proposed Smart-X mission would provide 30 
times the effecti ve area of Chandra w hile maintaining 
0.5" resolution (Vikhlinin et al.| 2012), making it the 
most ideal laboratory for probing X-ray scattering 
echoes. 
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